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ABSTRACT 

We present Spitzer 70/im and 160/im observations of the Spitzer extragalactic First Look Survey 
(xFLS). The data reduction techniques and the methods for producing co-added mosaics and source 
Qh catalogs are discussed. Currently, 26% of the 70/mi sample and 49% of the 160/im-selected sources have 

redshifts. The majority of sources with redshifts are star-forming galaxies at z < 0.5, while about 5% 
have infrared colors consistent with AGN. The observed infrared colors agree with the spectral energy 
distribution (SEDs) of local galaxies previously determined from IRAS and ISO data. The average 
160/im/70/im color temperature for the dust is Td — 30 ± 5K, and the average 70/im/24/im spectral 
index is a ~ 2.4 ± 0.4. The observed infrared to radio correlation varies with redshift as expected 
out to z ~ 1 based on the SEDs of local galaxies. The xFLS number counts at 70/im and 160/im are 
consistent within uncertainties with the models of galaxy evolution, but there are indications that the 
current models may require slight modifications. Deeper 70/im observations arc needed to constrain 
the models, and redshifts for the faint sources are required to measure the evolution of the infrared 
luminosity function. 

<^> Subject headings: galaxies: evolution — galaxies: starburst — infrared: galaxies 
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1. INTRODUCTION of depth or area coverage, the field has a large number 

^ r c , , ,. ■ j i • 1 1 i-i (~ 3000) of spectroscopic redshifts (P. Choi et al. 2005, 

One of the first science observations carried out with the > t-< i , , ■ , • \ 

5— U c m- c rp , /ttt , i r, nnA \ ,1 in preparation; h . Vlarleau et al. 2005, m preparation), 

i Spitzer Space lelescope (Werner et al. 2004) was the non- ,. //-■ , , , i ,• i • 
^ , . , , i i- t?- j. t i o / mo \ w ancillary radio (Condon et al. 2003 , and optical imag- 
es proprietary extragalactic _birst Look Survey xi'LS) which . . , J ,„ , , v , , , ., , , , ., , „. 
<__, i j , ? , • , -i « r j i i ■ i mg data (Fadda et al. 2004) that permit detailed multi- 
was designed to characterize the infrared sky at previously b , ,\ , , . ' , . . , . T , . . 

i a v ^n, td a c ■ ■ c. + a wavelength studies over a relatively large area. In this 

unexplored sensitivities. 1 he IR A S mission first uncovered to „ , i 

. , r ■ r j i • i • • iv. i i paper, we measure the source counts and use trie available 

ww the presence of infrared luminous galaxies m the local urn- ^ f. L , , . „ , r ; . ,. , . 

/at , , , mo a o vc ii ino7-i i , l redshifts to constrain the rest-frame spectral energy distri- 

verse (Neugebauer et al. 1984; boifer et al. 1987), and the . ,. ,„„„ s . . . ,. . r , 

ran ■ c j /en, j- i innn onno td j- v.- 4- i butions (SEDs) and derive the average infrared properties 

ISO infrared Elbaz et al. 1999, 2002; Rodighiero et al. P V T J ° , 

1 onncN j ju j u -u- 4 j -n- 4 i, lor the x_b LS 70/im and loOum populations. A cosmology 

rN , 2005) and ground-based submilhmeter and millimeter ob- f JJ — 70 km s^ 1 !^ c _1 fl —03 and ft — 7 is 

■ servations (Blain et al. 2002 and references therein) have a -s ° cd th o ho t tl ' ' e M ~~ ' ' A — • 

highlighted the importance that infrared luminous galax- assume roug ou l paper, 
ies have on the general understanding of galaxy evolution. 2 OBSERVATIONS 

Since the cosmic infrared background (CIB) peaks in the 

far-infrared (FIR) (Hauser & Dwek 2001), studying the The xFLS survey covers a 4deg 2 region (17 h 18 m 00 s , 

properties of galaxies that are bright in the FIR is crucial +59°30'00") within the northern continuous viewing zone 

for constraining models of galaxy evolution. of Spitzer 5 . Inside the xFLS main field a smaller verifica- 

In this paper, we present 70/im and 160/im observations tion strip of 0.25 deg 2 centered at 17 h 17 m 00 s , +59°45'00" 

of the xFLS field using the Multiband Imaging Photometer was observed with an integration time of 4 times that 

for Spitzer (MIPS, Rieke et al. 2004). The xFLS is a 4 deg 2 of the main survey to characterize the completeness and 

survey. The SWIRE Spitzer survey (Londsdale et al. 2004) s °urce reliability of the main survey. In total, 27.7 hours 

covers a wider area (49 deg 2 ) to similar depths, and deeper of xFL S observations were taken in 2003 December with 

observations covering smaller areas are being taken by the th e MIPS instrument. An additional 16.8 hours of ob- 

MIPS Instrument Team as part of the Guaranteed Time servations within the xFLS main field were taken in 2005 

Observers (GTO) program (e.g., Dole et al. 2004a) and May to characterize the performance of the 70/im array at 

other groups. Although the xFLS is not unique in terms warmer telescope temperatures (T mirror ~ 9.5 K). These 

observations produced useful data at 70/im, but the 160/im 

1 Spitzer Science Center, California Institute of Technology 220-06, Pasadena, CA 91125; frayerOipac. caltech.edu. 
2 Institut d' Astrophysique Spatiale, bat 121, Universite Paris Sud, F-91405 Orsay Cedex, France. 
3 Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721. 
4 Astronomy Program, Seoul National University, Shillim-dong, Kwanak-gu, Seoul, S. Korea, 2-880-9010. 
5 The extragalactic FLS data can be retrieved from the Spitzer Science Center at http://ssc.spitzer.caltech.edu/fls/ 
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data were not usable. Figure 1 shows the layout of the ob- 
servations. 

All of the MIPS observations were taken using the 
medium scan-rate mapping mode (4.2 s data collection 
events [DCEs]). The main-survey data were taken with 
adjacent 2deg scan legs that were offset in the cross-scan 
direction by 276" (nearly the full field of view of the ar- 
rays). The main-survey was covered twice to identify po- 
tential asteroids and to increase the redundancy of the 
data set. Unfortunately, half of the MIPS-70 array was 
rendered useless after launch due to a cable failure outside 
the instrument (Rieke et al. 2004). Therefore, at 70/im 
each position of the main survey was covered by only one 
scan leg. Since the data were not scheduled exactly as orig- 
inally planned, slight rotations between the 8 astronomical 
observation requests (AORs) yielded small gaps with zero 
coverage in the south-west corner of the main-survey 70/im 
map. 

The verification strip was observed using 4 AORs with 
0.5 deg scan legs at the medium scan rate. Each AOR 
covered the 0.5 deg x 0.5 deg field using cross-scan steps of 
148" (slightly less than half an array field of view). The 
warm test data taken in 2005 May were centered on the 
70/im main-field. These data consisted of 6 AORs with 
1.75 deg scan legs. Cross-scan steps of 148" were used to 
map the 1.75 degx 1.75 deg field once. Table 1 shows the 
average integration times, sensitivities, and area coverage 
for the 70/im and 160/im observations. 

3. DATA REDUCTION 
3.1. BCD Pipeline Processing 

The raw 70/im and 160/im (MIPS-Germanium [Ge]) 
data were downloaded from the Spitzer Science Center 
(SSC) archive and were reduced using the offline Ge Repro- 
cessing Tools (GeRT), available from the public SSC web- 
site. The GeRT uses an offline version of the SSC pipeline 
to produce the basic calibrated data products (BCDs), fol- 
lowing the algorithms derived by the MIPS Instrument 
Team and the MIPS Instrument Support Team (Gordon 
et al. 2005). The processing was done using the latest 
software (SSC pipeline version S12) to take advantage of 
improvements not currently available for the online xFLS 
data products (made with pipeline version Sll). 

BCD processing has two main steps: (1) calculation of 
the slope of the data ramp, and (2) calibration of the slope 
image. For the xFLS MIPS-Ge data, the raw 4.2 s data 
ramps are comprised of 32 non-destructive reads per pixel 
(one DCE). After correcting for the electronic nonlinearity, 
cosmic ray events and other discontinuities are identified 
in the data ramps. The SSC pipeline identifies disconti- 
nuities using a maximum likelihood technique (Hcssclroth 
et al. 2000). Linear slopes are calculated for the ramp 
segments between discontinuities and checked for consis- 
tency. The slopes from the segments are combined based 
on the empirical errors estimated from the scatter of the 
data within the ramp segments. Average slopes for each 
pixel are calculated to produce an uncalibrated slope im- 
age. 

The second step of BCD processing is the calibration of 
the slope image. The calibration of MIPS-Ge data is based 
on frequent measurements of internal stimulator flashes 

6 http://ssc. spitzer. caltech.edu/mips/dh/ 



(stims) which are used to track the responsivity of the de- 
tectors as a function of time. For the xFLS medium scan- 
rate data, the stim flashes are observed every 25 DCEs 
(118s). The stim flash signal is measured by subtract- 
ing the previous "background" DCE from the stim frame 
which is taken at the same position on the sky. For each 
AOR, a stim response function is calculated from 

interpolating between the stim minus background mea- 
surements. 

After the determination of the stim response as a func- 
tion of time (SR[t}), the BCD data are calibrated using 
the following equation: 

BCD(t) = FC[U(t)/SR(t) - DARK]/IC, (1) 

where U(t) is the uncalibrated slope image, DARK is the 
dark calibration file, and IC is the illumination correction 
calibration file which corrects for the combined illumina- 
tion pattern from the telescope and the stim flash signal. 
The DARK and IC calibration files are stable and are gen- 
erated by combining data from several different campaigns 
to improve the signal-to-noise. The flux conversion factor 
(FC) converts the instrument units into physical surface 
brightness units of MJysr -1 and is derived from observa- 
tions of standard calibrators (Sec. 3.4). 

3.2. BCD Filtering 

Before the 70/zm bias level was lowered in 2004 March, 
the 70/xm data showed significant data artifacts. Exam- 
ples and discussion of MIPS-Ge artifacts are shown in the 
MIPS Data Handbook 6 . The two main artifacts impact- 
ing the xFLS data are the stim flash latents and the vari- 
ations of the slow response as a function of time. The 
160/xm data are affected by these issues to a lesser degree 
due to the faster time constants of the 160/im stressed- 
Germanium detectors. For point sources, the stim latent 
and slow response residuals are additive effects. The fil- 
tering of the 160/im data is straight-forward. We used the 
filtered BCDs products, which remove a running median 
per pixel as a function of time by subtracting the median 
value for the surrounding DCEs closest in time (ignoring 
the current DCE, stim DCEs, and bad data). This high- 
pass time median filter removes data artifacts as well as 
the extended cirrus emission. To remove the stim latent 
artifacts at 160/im, we simply ignored the first DCE after 
the stim flash since the stim latents decay away within one 
DCE. 

The filtering process for the 70/im data is slightly more 
complicated. The time median filter by itself does not 
remove all of the data artifacts. Stim latents remain for 
many DCEs and are correlated by column. Since the scan 
map direction is nearly along the columns of the array 
(along the y-axis), the column artifacts are amplified. We 
remove the column residuals by subtracting the median of 
the values along each column for every BCD. The combi- 
nation of the high-pass time median filter per pixel and the 
column median filter removes the bulk of the data artifacts 
at 70/im. 

The resultant rms is lower for narrow high-pass median 
filter widths (< 15 DCEs), but narrow filter windows can 
yield significant negative side-lobes around bright sources. 
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In addition, column filtering does not maintain point 
source calibration for the brightest sources (j> 0.5 Jy). To 
avoid significant negative side-lobes and to preserve cal- 
ibration, the filtering of the MIPS-Ge data was done in 
two passes. The data from the first filtering pass were 
co-added, and sources were extracted (Sec. 3.3) to find 
the location of the bright sources. The source positions 
within the original BCDs were masked and new filtering 
corrections were calculated in the second pass, ignoring 
the pixels containing sources. This second-pass filtering 
technique minimizes the data artifacts while preserving 
point-source calibration. 

Different filtering methods were tested and optimized 
using the GeRT. At 70^m, we applied column filtering fol- 
lowed by a high-pass time median filter with a width of 
12 DCEs to yield the deepest image. The final sensitivity 
of the filtered 70/im image was improved by more than 
a factor of two in comparison to data with no filtering. 
Filtering has less of an impact at 160/xm, but is useful 
for deep observations of faint point sources. We adopted 
a high-pass time median filter width of 12 DCEs for the 
160^m data using the second-pass filtering technique as 
done for 70/xm. 

3.3. Data Coaddition and Source Extraction 

The filtered BCDs were coadded using the SSC mosaic- 
ing and source extraction software (MOPEX) (Makovoz 
& Marleau 2005). The data were combined ignoring bad 
data flagged in the BCD mask files (bmask) and the bad 
pixels defined by the static pixel mask (pmask). Using the 
redundancy of the MIPS data, additional spurious data 
values were removed via outlier rejection. The remaining 
data were corrected for array distortions and projected 
onto a sky grid with square pixels that were more than 
2 times smaller than the original pixels. The output mo- 
saic pixels are 4" at 70^m and 8" at 160^m. The data 
were averaged using weights proportional to the fractional 
area subtended by the original pixels as projected onto the 
output pixel grid. Since the calibration of the main, verifi- 
cation, and the warm 70/zm data sets are consistent within 
measurement errors, all data were combined to yield the 
deepest images. 

Sources were extracted from the final images using 
MOPEX. The MOPEX software uses a Point-source Re- 
sponse Function (PRF) image to optimize source detec- 
tion and point source fitting. A PRF calibration image 
was made for each MIPS band by coadding isolated bright 
point sources within the xFLS field. A point source prob- 
ability (PSP) image was made from non-linear matched 
filtering of the input image with the PRF. The PSP image 
represents the probability of having a point source above 
the noise at each pixel. The PSP image was multiplied by 
the input image to yield the filtered image used for source 
detection. The filtered image enhances the presence of 
point sources in the mosaic while smoothing out noise fea- 
tures that do not match the input PRF, providing a more 
robust image for source detection. After source detection, 
the original images were fitted with the PRF image via x 2 
minimization to extract source positions and flux densi- 
ties. Regions containing multiple peaks were de-blended. 
For extended sources and/or close blends, not well fitted 

7 The effective wavelengths of the MIPS-24, 70, and 160 bands are 



by the PRF, large aperture measurements were used to 
derive the total flux density for the source. 

The noise across the images is non-uniform. For optimal 
source detection and extraction, an accurate representa- 
tion of the uncertainties is needed. The MOPEX uncer- 
tainty image based on the input BCD pipeline uncertain- 
ties is useful for identifying the relative uncertainties across 
the image, but the pipeline errors typically underestimate 
the absolute uncertainty level for low background regions. 
The median level of the MOPEX uncertainty image was 
scaled to match the average empirical noise in the mosaics. 
By using a properly scaled uncertainty image, we find an 
average x 2 f° r the fitted point sources of x 2 ~ 1, suggest- 
ing that the uncertainty image is self-consistent with the 
average point source errors in the data. 

Several tests were done to validate the photometric re- 
sults from MOPEX. No calibration differences were found 
for different sub-sets of the data. We also compared the 
photometry done with the MOPEX software with that 
from the program StarFinder (Diolaiti et al. 2000). Both 
programs use an empirical PRF to fit source flux density. 
The results from MOPEX and StarFinder agree to better 
than 5% for bright point sources (;> 50mJy at 70^m and 
^>200mJy at 160^m) and better than 10% rms down to 
the limits of the catalog, consistent with the uncertainties 
in the data. For faint sources detected in the MIPS-Ge 
bands, PRF fitting techniques yield more reliable results 
than aperture measurements. By comparing results be- 
tween the verification data and the main survey, we find a 
consistency of ±10% for PRF matching techniques and a 
dispersion for aperture measurements of larger than ±25% 
at faint flux density levels. 

3.4. Calibration 

The absolute calibration is derived using observations 
of stars at 70^m and observations of asteroids and well- 
studied luminous infrared galaxies at 160^im. The data 
were processed assuming flux conversion factors of 634 
MJysr" 1 per MIPS-70 data unit and 42 MJysr" 1 per 
MIPS-160 data unit. The flux densities in the source cat- 
alogs were multiplied by the correction factors of 1.16 and 
1.17 for the MIPS-70 and MIPS-160 bands respectively, 
which includes the color corrections associated with a con- 
stant vj v SED and updated calibration conversion factors 
based on the latest measurements. Additional color cor- 
rections for sources with non-constant vj v SEDs are negli- 
gible since the color corrections are similar for the range of 
SEDs appropriate for galaxies; f v oc v~ a , where a ~ 0- 
3. For flux density ranges of 50mJy-2Jy, the absolute 
calibration uncertainty is estimated to be about 15% and 
25% for the 70/xm and 160/im bands respectively. No cor- 
rections have been made for possible flux nonlincarities, 
which may be significant for sources brighter than this 
nominal flux density range. Observers are recommended 
to check the latest information on MIPS calibration at the 
SSC web site. 

We confirmed that the absolute 70/im flux density scale 
for the xFLS data is consistent with previous measure- 
ments made by IRAS to within 10%. For comparison with 
IRAS, we interpolated between the flux densities in IRAS- 
60 and IRAS-100 bands (Beichman et al. 1988) to the ef- 

■.7/im, 71.4/im, and 155.9/^m, respectively. Throughout this paper, the 
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fective wavelength of the MIPS-70 band (71.4/im) 7 apply- 
ing the appropriate color corrections. For the five bright- 
est sources in the field detected by IRAS, we derive a flux 
density ratio of Spitzer/IRAS = 1.0 ± 0.1 at 71.4/im. The 
calibration at 160/im is more uncertain given that xFLS 
sources have not been observed previously at wavelengths 
longer than 100/im. Simple SED model fits to the data 
indicate that the calibration at 160/im is consistent with 
IRAS to within 30%. 

3.5. Mosaics and Catalogs 

The MIPS-Ge mosaics and catalogs are available online 
at the SSC website. For both of the MIPS-70 and 160/mi 
bands, we provide coadded mosaics of the entire xFLS data 
sets along with the corresponding coverage maps and un- 
certainty images. The mosaics are from the filtered prod- 
ucts and have been background subtracted. The absolute 
level of the uncertainty image has been scaled to match the 
empirical noise within the mosaics. The coverage map rep- 
resents the effective number of input BCDs for each pixel 
of the mosaic after outlier rejection. There are regions 
within the 70/im and 160/im mosaics with zero coverage. 
The 70/xm image has gaps in the south-west corner of the 
mosaic due to non-overlapping AORs. The 160/im mosaics 
have a regular pattern of low coverage due to the unusable 
block of five detectors and the low redundancy of these 
observations. 

Point-source catalogs were made for sources with S/N> 
7 at each wavelength to insure high reliability. Sources 
near the edges or within regions of very low coverage were 
deleted from the public catalogs to avoid potentially spu- 
rious sources. In total, 687 70/zm sources (Table 2) and 
207 160/xm sources (Table 3) are cataloged. The 70/xm 
and 160/im catalogs are independent. The catalogs can 
be combined together or with xFLS data at other wave- 
lengths, depending on the specific application. The aver- 
age point source rms (lex) is 2.8 mJy and 20mJy for the 
70/im warm+main data and 160/im main-survey data, re- 
spectively. For the deeper verification regions, the point 
source rms is 1.6 mJy at 70/xm and lOmJy at 160/im. 

The images were examined visually to validate source 
detection and to determine the appropriate method for 
deriving the flux densities on a source by source basis. 
About 98% of the derived flux densities are based on PRF 
fitting (Sec. 3.3). The remaining sources mostly consist 
of bright extended galaxies and/or close blends that are 
not well fitted by the PRF. In these cases, large aperture 
measurements were made to derive the total flux densities. 
For bright extended sources that have one or more neigh- 
boring point sources within the aperture, the flux density 
of the bright galaxy was derived by subtracting the PRF 
measurement (s) of the faint point source(s) from the total 
aperture measurement. 

Table 2&3 show the format for an example portion of 
the xFLS MIPS 70/im and 160/zm catalogs published in the 
online edition of the Journal. The average radial positional 
errors are 2'.'6 and 5'/ 2 (la) for the S/N> 7 MIPS-70 and 
160/im sources respectively. No systematic differences are 
found for the 70/im source coordinates (< 0'.'2) in compar- 
ison to the more accurate 24/im positions. For consistency 
with the 24/zm and 70/zm data, we correct the 160/im po- 

flux densities for the MlPS-bands are defined as S24= S„(23.7/tm), S70= 



sitions by 4'.'7 to compensate for the systematic positional 
offset measured for the 160/zm sources (a known issue for 
pre-S12 versions of the pointing pipeline). The flux den- 
sities have been color corrected assuming a constant vf v 
SED. The tabulated flux density errors include the abso- 
lute flux density uncertainties of 15% and 25% at 70/im 
and 160/im, respectively. 

3.6. Reliability and Completeness 

The verification data were used to test the reliability of 
the source detection and extraction techniques. For a high 
level of reliability at S/N < 10, it is important that the 
uncertainty image accurately reflects the small-scale spa- 
tial variations in the noise across the images. Using the 
uncertainties based on the BCDs, properly scaled to rep- 
resent the average empirical noise in the data (Sec. 3.3), 
we obtained good results. Based on the deeper verifica- 
tion data, we do not find spurious sources with S/N > 5 
in the main survey within the verification field. Sources 
detected at S/N > 4 are 80-85% reliable, and detections at 
the S/N > 3 level are only 50-60% reliable. Based on these 
results, we adopted a S/N > 5 cut for deriving the source 
counts in Sec. 4.1. A more conservative S/N> 7 criterion 
was adopted for the public catalogs (Sec. 3.5). 

At 70/im the completeness and source counts were mea- 
sured within the overlapping region of the warm and main 
fields (Fig. 1), while the entire main-field xFLS area was 
used at 160/im. Both empirical completeness estimates 
based on the verification data and simulated completeness 
measurements were made. The simulations of complete- 
ness were carried out by adding point sources with differ- 
ent flux densities into the mosaics at random locations, and 
then extracting sources using the same techniques adopted 
for the source catalogs. The simulations and empirical 
methods for estimating completeness gave consistent re- 
sults. Using the adopted S/N> 5 criterion, the complete- 
ness level falls rapidly below 60% for S70< 14mJy (Fig. 
2) and S160< 100 mJy (Fig. 3). Similar completeness sim- 
ulations were done for the deeper verification field, and we 
find 60% completeness levels of approximately 9 mJy and 
60mJy at 70/im and 160/im, respectively. 

4. RESULTS AND DISCUSSION 
4.1. Source Counts 

The source counts are derived for the main xFLS and 
verification fields based on S/N> 5 catalogs corrected for 
completeness. At 70/im we use the 3.3 deg 2 region con- 
taining both the main and warm survey data. The 160/im 
counts are based on the entire 4.5 deg 2 area of the main 
survey. Figures 4&5 show the differential source counts 
(dN/dS x S 2 - 5 ) at 70/im and 160/im respectively. Over 
the entire observed range of flux density, both the 70/xm 
and 160/xm counts increase at super-Euclidean rates with 
decreasing flux density. 

In total 845 sources at 70/im were detected (> 5cr at 
S70< 440 mJy) within the main+warm field area and 186 
sources were detected (S70< 70mJy) within the verifica- 
tion field (Table 4). In comparison at 160/im, 227 sources 
were detected in the main-field (S160< 880 mJy) and 45 
sources (S160< 140 mJy) were found in the verification 
field (Table 5). The source counts are consistent between 

(71.4/im), and S160= S„(155.9/im). 
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the main and verification data for the overlapping range of 
flux densities. The xFLS counts are also consistent within 
errors with the previous MIPS-Ge counts published by 
Dole et al. (2004a). The error bars for the source counts 
include the Poisson noise, the completeness errors, and the 
absolute flux density calibration uncertainty. The data 
were binned such that the errors are not dominated by 
Poisson statistics, except for the highest flux density bins. 
In the lowest flux density bins, the errors on completeness 
dominate the uncertainties. 

Dole et al. (2004a) found a difference in the 70^m counts 
between the CDF-S and Marano observations, which they 
attributed to field variations. The measured 70^m xFLS 
counts are within the range of values previously measured, 
and closer to the CDF-S counts at faint flux densities. The 
xFLS-Ge counts are also consistent within uncertainties 
with the evolutionary models of Lagache et al. (2004), 
which were updated to match the Spitzer 24,um counts. 
However, the xFLS data points are systematically slightly 
lower than the predictions from current models. Counts 
from the wide-area SWIRE survey would be needed to 
check for possible field variations and to test whether the 
models require slight modifications. 

Although the xFLS 70/xm counts are at fainter levels 
than those previously published (Dole et al. 2004a), they 
are not yet deep enough to measure the location of the 
expected turn over of the differential counts (Fig. 4) . The 
flux density at which the differential counts turn over pro- 
vides important constraints on evolutionary models (La- 
gache et al. 2004; Chary & Elbaz 2001; King & Rowan- 
Robinson 2003; Xu et al. 2001), as shown by the 24^m 
data (Marleau et al. 2004; Papovich et al. 2004; Chary 
et al. 2004). At 160/im Spitzer is not expected to mea- 
sure the turn over in the counts (Fig. 5), since the 160^m 
data are expected to be limited by confusion below 40 mJy 
(Dole et al. 2004b). The differential counts at 70^m are 
expected to turn over at around lOmJy (Lagache et al. 
2004), which is well above the confusion limit, suggesting 
that deeper 70^m counts would provide useful constraints 
for galaxy evolution models. 

4.2. Infrared Colors 

For high reliability, the source counts were derived using 
a S/N> 5 criterion for the individual bands. However, it 
is possible to maintain a high level of reliability at lower 
S/N by using the MIPS-24/im and/or radio data of the 
field. For studying the infrared colors, we used S/N> 4 
source lists to increase the number of sources detected at 
high-redshift in the MIPS-Ge bands. The 70/xm source 
list was bandmerged with the 24/xm catalog (D. Fadda et 
al., in preparation). At S/N> 4, 1779 xFLS-70 sources 
have counterparts at 24^m within 6". The different reso- 
lutions of the 24^m, 70/xm, and 160/^m data sets (6", 18'.' 5, 
and 40", respectively) complicate the identification of the 
counterparts in the MIPS bands. Out of the 1779 70/im 
sources, 9% have multiple candidate 24/iui counterparts. 
To avoid potential spurious matches, we only used sources 
with one-to-one matches for studying the infrared colors. 
The 4er xFLS-160 source list was matched to the 70^m 
data set using a positional matching radius of 13". The 
resulting matched MIPS-160+70 source list was matched 

8 We use 7 to represent the "a" parameter of Dale et al. (2001). 



to the MIPS-70+24 source list. We find 301 160^m sources 
with one-to-one matches between all three MlPS-bands. 

Out of the 1618 xFLS-70 sources with one 24/um posi- 
tional match, 427 currently have redshifts. Approximately 
57% of these sources come from the NOAO WIYN-Hydra 
radio-selected redshift survey (F. Marleau et al. in prepa- 
ration), 16% of the redshifts are from the Keck-DEIMOS 
24^m-selected sample (P. Choi et al., in preparation), and 
the remaining 27% of redshifts are from the SLOAN survey 
(Strauss et al. 2002). The largest redshift to date is for a 
quasar at z = 3.56. For the 160/zm-selected sources with 
detections in all three MlPS-bands, about half (146) have 
known redshifts. Although the redshift survey of xFLS- 
Ge sources is not nearly complete (only 26% of the 70^m 
sources and 49% of the 160/xm sources currently have red- 
shifts), there are a sufficient number of sources to study 
the general trends in the infrared colors as a function of 
redshift for the sample. 

We adopt two sets of SED models for studying the in- 
frared colors. The first set of SED models is based on 
a simple modified blackbody (Blain, Barnard, & Chap- 
man 2003). The SED is expressed as /„ — e u B u , where 
B v {Td , v) is the blackbody function for a dust temperature 
Td, and e„ oc v" is the dust emissivity. In the mid-IR, we 
substitute a power-law of the form /„ cx v~ a , smoothly 
matching e v B v at longer wavelengths (Blain et al. 2003). 
For this simple SED model, there are three free param- 
eters: Td, P, and a. The S160/S70 ratio constrains Td, 
while the S70/S24 ratio measures a. The dust emissivity 
index (3 is not very well constrained by the Spitzer data 
alone. For simplicity, we adopt a constant value of (3 = 1.5 
which is consistent with the results for low-redshift IRAS 
galaxies (Dunne et al. 2000). 

The second set of SED models uses a physically more 
realistic approach by assuming a power-law distribution of 
dust masses as a function of radiation field intensity (Dale 
et al. 2001). The dust mass {Md) as a function of the 
radiation field (U) is given by dMd{U) oc U~ 7 dU (Dale et 
al. 2001; Dale & Helou 2002) 8 . Quiescent, cirrus-like dust 
regions are expected to have 7 ~ 2.5, while environments 
near active Hn regions are expected to be approximated 
by 7 ~ 1. For a mixture of active and quiescent regions, 
the average effective 7 of star-forming galaxies should be 
between 1 7 <J 2.5. 

Figure 6 shows the observed 70/zm/24^m flux density ra- 
tio (S70/S24) plotted as a function of redshift with several 
predictions based on the SEDs of local galaxies. The ma- 
jority of the xFLS sources have S70/S24 ratios within the 
range of values between the extreme ultraluminous star- 
burst of Arp200 and the more typical starburst M82. At 
moderate and high redshifts, infrared-cool ultraluminous 
infrared galaxies (ULIRGs, L > 10 12 L Q ), e.g., Arp220, 
are expected to have the highest S70/S24 ratios, while 
warm-infrared ULIRGs, e.g., Mrk 231, and AGN domi- 
nated sources are expected to have lower S70/S24 ratios. 
Spiral galaxies without a strong infrared excess in their 
SED (e.g., M100) show decreasing S70/S24 ratios as a 
function of redshift; however, these galaxies are not de- 
tected at high redshifts (z > 0.5) in the xFLS-Ge survey. 
In general, the S70/S24 ratio can be used to help distin- 
guish between AGN-dominated sources from star-forming 
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galaxies, except in cases for which strong poly cyclic aro- 
matic hydrocarbon (PAH) emission features are redshiftcd 
into the 24/im band (e.g., M82 at z ~ 2 in Fig. 6) and for 
AGN with strong silicate absorption features (e.g., Spoon 
et al. 2004) at z ~ 1.5. 

For the xFLS sample currently with redshifts, only 5% 
show low S70/S24 ratios consistent with AGN SEDs. The 
AGN population is fitted with T d = 90 ± 30 K and a = 
1.1 ± 0.3 (Fig. 6), which is consistent with the IRAS ob- 
servations of Seyfert galaxies (Miley, Neugebauer, & Soifcr 
1985). The majority of the xFLS sample with redshifts are 
starbursts at z < 0.5. Figure 6 shows a large concentration 
of sources at z <~ 0.2. The S70/S24 ratios based on the 
SEDs of quiescent and active galaxies (Dale & Helou 2002) 
intersect at z ~ 0.2, which may contribute to this concen- 
tration of data points. For the galaxies with z ~ 0.2, the 
average observed infrared color is S70/S24= 14 ± 5. In 
comparison with the IRAS population (Soifer et al. 1989; 
Sanders et al. 2003), the IRAS-60 / IRAS-25 ratios corre- 
spond to a predicted ratio of S70/S24= 14 ± 3, consistent 
with the measured xFLS infrared colors. For the simple 
SED model, the average S70/S24 ratio corresponds to an 
infrared index of a = 2.4. Using SEDs of Dale & Helou 
(2002), values of 7 ~ 2 give the best agreement with the 
average S70/S24 ratios. 

Although the majority of the 70/xm sources do not have 
redshifts, the distribution of the S70/S24 ratios is the same 
for galaxies with and without redshifts (Fig. 7). The av- 
erage mid-IR spectral index is a = 2.4 ±0.4. For a simple 
power-law representation of the mid-IR SED, the modeled 
S70/S24 ratio is constant with rcdshift (e.g., Starburst 
solid line in Fig. 6). However, galaxies are expected to 
show significant variations in the S70/S24 ratio as strong 
PAH emission and mid-IR absorption features (e.g., Ar- 
mus et al. 2004) are redshifted through the 24/um band 
(as seen for Arp 220 in Fig. 6). 

Using simple modified blackbody SEDs, we calculate 
S160/S70 as a function of redshift and dust temperature 
for a constant value of a = 2.4 (Fig. 8). The average 
160^m/70yum color temperature for the dust temperature 
is 30±5K. This temperature is slightly lower than the av- 
erage temperature of T d — 38 ± 3K (Dunne et al. 2000) 
and T d ~ 35-40 K (Soifer et al. 1989) derived for IRAS 
galaxies, but the longer wavelength 160^m band is more 
sensitive to regions with cooler dust temperatures than 
those observed by IRAS. The blackbody temperature of 
30 K derived here agrees well with the long- wavelength ISO 
observations of infrared-bright spiral galaxies (Bendo ct al. 
2003). Bendo et al. (2003) found that T d ~ 30 K provided 
the best fit to the ISO data for (3 ~l-2, consistent with 
the results in Figure 8. In the context of the SED models 
of Dale & Helou (2002), values of 7 ~ 2-2.5 provide the 
best agreement with the observed S160/S70 ratios. 

The distribution of the S160/S70 ratios is measurably 
different for sources with and without redshifts (Fig. 9). 
Unlike the S70/S24 ratio which is on average roughly con- 
stant with redshift, the S160/S70 ratio is expected to in- 
crease with redshift (Fig. 8). Figure 9 shows an excess 
of 160/zm sources without redshifts having high S160/S70 
ratios (Log[S160/S70]^ 0.8). These results may suggest 
that the excess sources are at high redshift (z J> 0.5) and 
that the current redshift sample of 160/iin sources may be 



biased toward lower redshifts. 

4.3. Infrared Luminosity 

For the 70^m selected sources, we can estimate the av- 
erage infrared luminosity for the xFLS population. The 
rest-frame wavelength of the MIPS-70 band corresponds to 
60/im at a redshift of z ~ 0.2. The FIR luminosity (42.5- 
122.5/um) can be derived from the rest-frame 60^m and 
100/im flux densities (Helou, Soifer, & Rowan-Robinson 
1985). By using the SED associated with the average 
model parameters of T d ~ 30 K and a ~ 2.4 for the xFLS 
sources, we estimate a rest-frame S100/S60 ratio of 2.3. 
This is consistent with the values of S100/S60= 2 ± 1 ob- 
served for the bulk of the IRAS galaxies (Sanders et al. 
2003). Based on the rest-frame S100/S60 ratio, the corre- 
sponding FIR flux is FIR(Wm- 2 ) = 6.1 x 10~ 14 (S70/Jy) 
for the observed S70 flux density at z ~ 0.2. For galaxies 
at z ~ 0.2, the average value of S70 is 33 mJy This cor- 
responds to a FIR luminosity of L(FIR)= 6.0 x 10 10 L & . 
Using the models of Dale et al. (2001) and the estimated 
S100/S60 ratio, the total infrared luminosity (3-1100/um) 
L(TIR)= 2.3L(FIR). Hence, the average total infrared lu- 
minosity of the 70/xm-selected xFLS galaxies at z ~ 0.2 is 
about L(TIR)= 1.4 x 10 11 L Q . 

Although the redshift survey for the 70yum-selected 
sources is not currently complete, the infrared luminosity 
function can be estimated at low redshift for bright flux 
densities. At S70> 50mJy, 72% (65/90) of sources have 
redshifts. The majority of sources without redshifts are 
suspected to be at high redshift given that the observed 
rcdshift distribution declines significantly at z > 0.3 and 
the models predict that the majority of sources at these 
flux density levels are at z > 0.3. Based on the mod- 
els of Lagache et al. (2004), about 40-45% of galaxies 
with S70> 50mJy are predicted to be at z < 0.3. As- 
suming ~ 100% completeness at z < 0.3 (i.e., assume 
all sources with S70> 50mJy currently without redshifts 
are at z > 0.3), we find 64±9% (58/90) of galaxies with 
S70> 50mJy and z < 0.3. If the current xFLS redshift 
surveys arc not complete at z < 0.3 (with S70> 50mJy), 
then the models would be even more discrepant with the 
observations. Hence, the models may under predict the 
percentage of low redshift 70/xm sources at bright flux den- 
sities. There are two potential caveats on this result. First, 
a small percentage (< 1%) of 70/xm sources do not have 
24/im counterparts, and these sources are likely to be at 
high redshift (e.g., at z ~ 1.5 where the silicate absorption 
feature is redshifted into the 24^m band). However, there 
is only one 70/xm source without a 24/xm counterpart with 
S70> 50mJy, so this population, by itself, cannot account 
for the apparent discrepancy with the models. The sec- 
ond caveat is that we have thrown out 9% of the 70/xm 
sources with multiple possible 24yum matches (Sec. 4.2). 
The 70/im sources with multiple candidate 24/im counter- 
parts tend to be brighter on average and typically have low 
redshifts. In fact, we find only one such source with z > 0.3 
and with S70> 50mJy (1/70 with redshifts). Additional 
rcdshift surveys are needed to confirm the findings here 
which suggest that the current models slightly over pre- 
dict the percentage of bright 70/zm sources (S70> 50 mJy) 
at high-redshift (z > 0.3) and their contribution to the 
total CIB. 

To estimate the infrared luminosity function, we use the 
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sample of 58 galaxies with S70> 50mJy and z < 0.3. 
We adopt the 1/V max method following the calculation 
of the IRAS vLv(60pm) luminosity function (Soifer et al. 
1987). For each source, the maximum volume (V max ) at 
which the source would be included in the sample is com- 
puted, using V max (z — 0.3) as an upper bound for the sur- 
vey. The 70/im flux density is converted to a rest-frame 
60/im flux density using a power-law interpolation between 
the observed 24/xm and 70/im measurements. The calcu- 
lated xFLS 60/im luminosity function (p) shown in Fig- 
ure 10 represents the space density of galaxies per Mpc 3 
per magnitude in luminosity. Even though the current 
xFLS sample is a factor of 100 deeper in flux density 
than the IRAS 60/im bright galaxy sample (Soifer et al. 
1987; Sanders et al. 2003), the xFLS sample does not 
yet probe the high luminosity end of the luminosity func- 
tion (Log(i/L„ [60/im]/ Lq) > 11.5), given the low redshift 
range (median redshift of only z = 0.096) and small area 
of the current sample. The derived xFLS luminosity func- 
tion agrees well with the IRAS 60/im luminosity function 
(Soifer et al. 1987). Below Log(^L„ [60/im]/ Lq) < 10.5, 
p oc {y L„)~ - 8 (Soifer et al. 1987), while at higher lumi- 
nosities p oc (yL v )~ 2 - 2 (Sanders et al. 2003). Additional 
rcdshifts for the faint S70~10mJy sources are needed to 
constrain the evolution of the luminosity function out to 
z ~ 1. 

4.4. Infrared to Radio Correlation 

In the local universe, galaxies over a wide range of lu- 
minosity and Hubble types follow the empirical IR/radio 
correlation (Helou et al. 1985; Condon 1992). Observa- 
tions of ISO sources (Gruppioni et al. 2003) indicate that 
the mid-IR to radio relationship holds out to z ~ 0.6, 
and the data for the SCUBA sources suggest that the FIR 
to radio relationship may be applicable even at z ~ 2-3 
(e.g., Chapman et al. 2005). With Spitzerwe can measure 
the IR/radio relationship out to z ~ 1 by direct observa- 
tions near the peak of the SED, instead of estimating the 
IR luminosity from the shorter mid-IR or longer sub-mm 
wavelengths. Appleton et al. (2004) presented the first es- 
timate of the IR/radio relationship based on Spitzer data. 
Their results were based on the early analysis of the xFLS 
data. Here we present updated results using more sensitive 
70/im data. 

The sample of 427 xFLS-70 sources with redshifts and 
with only one candidate 24/im counterpart (Sec. 4.2) was 
matched to the 4<r radio catalog (Condon et al. 2003). 
Within 3", 325 of these xFLS-70 sources (76%) have radio 
counterparts. Figure 11 shows the observed S70/S(20cm) 
ratio as a function of redshift for the matched galaxies. 

At low redshift, the FIR to radio q-parameter is 
defined as q = Log{FIR/(3.75 x 10 12 Wm~ 2 )) - 
Log(S'(20cm)/(Wm" 2 Hz" 1 )) = 2.3 ± 0.2 (Helou et al. 
1985). Based on the typical SED of the xFLS sam- 
ple of galaxies (Sec. 4.2), the corresponding predicted q 
parameter for the 70/im to radio flux density ratio is 
<77o = Log(570/S'[20cm]) = 2.09 at z ~ 0.2. We measure 
an average value of q^o = 2.10 ± 0.16 for xFLS galaxies at 
z ~ 0.2. By comparison, Appleton et al. (2004) derived a 
similar value of q^o = 2.16 ± 0.17, but unlike the previous 
results we observe the expected decrease of observed (770 
as a function of redshift (Fig. 11). Based on the SEDs of 
local galaxies, we expect the observed 570 parameter to de- 



crease with redshift since the average IR spectral index is 
significantly steeper than the radio spectral index. Assum- 
ing an IR spectral index of a = 2.4 based on the observed 
average ratio of S70/S24= 14 and the typical non-thermal 
radio spectral index of a = 0.8 (/„ oc v~ a , Condon et al. 
1992), we expect qro oc (1 + z) -16 . From a least-squares fit 
to the data, we estimate q?o oc (1 + z)~ 1A±(l6 . Hence, the 
observed q?o for the xFLS galaxies follows the expected 
trend out to at least z ~ 1 (Fig. 11). 

5. CONCLUSIONS 

The Spitzer 70/im and 160pm observations of the xFLS 
are presented. For the deeper verification field data, we 
measure number counts down to 8mJy and 50mJy (5cr) 
in the 70/im and 1 60/im bands, respectively. The observed 
xFLS counts are consistent with previous measurements 
(Dole et al. 2004a) and are consistent within the un- 
certainties with the evolutionary models (Lagache et al. 
2004). Based on the models of Lagache et al. (2004), 
approximately 35% of the CIB is resolved at 70/im and 
15% is resolved at 160/im at the depth of the verifica- 
tion data. The observed fraction of low redshift galaxies 
at bright 70/im flux densities is larger than model predic- 
tions and the total counts appear slightly systematically 
lower. These results may suggest the models overestimate 
the contribution of high-redshift galaxies at bright 70/im 
flux densities to the total CIB. Deeper 70/im observations 
are needed for measuring the expected turn over in the 
differential source counts to provide better constraints on 
the evolutionary models. 

The observed xFLS infrared colors S70/S24 and 
S160/S70 are consistent with the results from the IRAS 
and ISO missions. Modeled SED fits suggest an average 
160/im/70/im color temperature for the dust of Td ~ 30 K, 
for the 160/im-selected sample of galaxies. This tempera- 
ture is consistent with ISO observations of spirals (Bendo 
et al. 2003). The average S70/S24 ratio implies an infrared 
spectral index of a ~ 2.4, which agrees with expectations 
from the average IRAS S60/S24 ratio (Soifer et al. 1989; 
Sanders et al. 2003). The observed 70/im infrared to ra- 
dio correlation of the xFLS sources also agrees well with 
the FIR-to-radio correlation found for local star-forming 
galaxies (Helou et al. 1985; Condon et al. 1992). We ob- 
serve a trend of decreasing S70/S(20cm) as a function of 
redshift consistent with expectations based on the SEDs 
of local galaxies. 

Given the lack of sufficient redshift measurements for 
the high-redshift faint xFLS-Ge sources, we can only de- 
rive an infrared luminosity function at low redshift for only 
the brightest 70/im galaxies. With these data, we calculate 
a rest-frame 60/im luminosity function for the xFLS that 
agrees well with the IRAS luminosity function (Soifer et 
al. 1987). In the future, it may be possible to measure the 
evolution of the luminosity function out to moderate red- 
shifts with the xFLS data when additional redshift surveys 
of faint 70/im sources become available. 

We thank all of our colleagues associated with the 
Spitzer mission who have made these observations pos- 
sible. This work is based on observations made with 
the Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology 
under NASA contract 1407. 
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Fig. 1. — The layout of the 70/tm and 160/im fields for the xFLS observations. The 70/im main and verification fields are shown by the 
solid-lines, while the 160/im xFLS fields are shown as dashed-lines. The location of the 70/im-warm observations taken in 2005 May is shown 
by the dotted-line. The right ascension and declination offsets are with respect to the center of the main field (17 h 18 m 00 s , +59°30'00"). 
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Fig. 2. — The 70/xm completeness level as a function of flux density for the verification and warm+main xFLS surveys. The diamonds 
represent empirical measurements of the warm+main survey based on the deeper verification data, and the squares and solid circles are the 
results of simulations for the warm+main and verification fields respectively. The completeness estimates arc for the S/N> 5 source lists used 
to derive the counts. 
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Fig. 3. — The WOfim completeness level as a function of flux density for the verification and main xFLS surveys. The diamonds represent 
empirical measurements of the main survey based on the deeper verification data, and the squares and solid circles are the results of simulations 
for the main and verification fields respectively. The completeness estimates are for the S/N> 5 source lists used to derive the counts. 
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Fig. 4. — The 70^im differential number counts for the verification (solid circles) and warm+main (squares) xFLS surveys. The diamonds 
represent previous MIPS-GTO measurements (Dole et al. 2004a). The vertical error bars include the uncertainties of the absolute flux density 
scale and completeness, and the horizontal line segments show the bin sizes. The solid line is the evolutionary model of Lagache et al. (2004). 
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Fig. 6. — The S70/S24 flux density ratio for 70^m sources with redshifts. The ULIRGs Arp220 and Mrk231, the starburst M82, and 
the spiral galaxy M100 are shown as dashed lines. The solid lines represent simple modified blackbody SED models for starbursts (a = 2.4, 
T d = 30 K) and AGN (a = 1.1, T d = 90 K). The dotted lines show the SED models of Dale & Helou (2002) as a function of 7, where 7 ~ 2-2.5 
represent more quiescent star-forming galaxies and 7 <~ 1.5—1 represent more active galaxies. The plus symbol shows the predicted S70/S24 
ratio at a rcdshift z = 0.2 based on the average IRAS S60/S25 ratio. The typical error bar for the data points from the 4cr source lists is 
shown in the lower right. 
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Fig. 7. — The distributions of the S70/S24 flux density ratios for 70/im sources with redshifts shown in Fig. 6 (dashed line) and for all 
70/im-sclcctcd sources (solid line). Values for a range of infrared spectral indexes (a) are shown by the vertical dotted lines. 
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Fig. 8. — The S160/S70 flux density ratio for 160/^m-selected sources with redshifts. The solid lines show simple modified blackbody SED 
models based on different dust temperatures (20-50 K) assuming an infrared spectral index of a = 2.4 and a dust emissivity index of (i = 1.5. 
The dotted lines show the SED models of Dale & Hclou (2002) as a function of 7. The typical error bar for the data points from the 4a 
source lists is shown in the lower right. 
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Fig. 9. — The distributions of the S160/S70 flux density ratios for the 160/im sources with redshifts shown in Fig. 8 (dashed line) and for 
all 160/im-selected sources (solid line). The excess of sources without redshifts at the higher S160/S70 ratios may suggest that the current 
redshift sample of 160/^m sources is biased toward lower redshifts. 
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Fig. 10. — The xFLS rest-frame 60^im luminosity function based on the sources with S70> 50mjy and z < 0.3 (solid circles). The IRAS 
60/im luminosity function is represented by the squares (Soifcr et al. 1987), corrected for the adopted value of Ho = 70km s — 1 Mpc -1 . The 
dottcd-linc shows the double power-law fit of the luminosity function based on the IRAS data (Soifer et al. 1987; Sanders et al. 2003). 
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Fig. 11. — The S70/S(20cm) flux density ratio for the 70/im-selected xFLS sources with redshifts between < z < 1. The dotted-line shows 
the expected behavior of S70/S(20 cm)oc (1 + 2)~ 1,6 , assuming power-law approximations for the SEDs at radio and infrared wavelengths. 
The dashed-line shows the result of a least-squares fit [S70/S(20cm)oc (1 + z)- 14±0 - 6 ]. The plus symbol shows the predicted value of the 
S70/S(20cm) ratio based on the local infrared to radio correlation. The typical error bar for the data points from the 4<r source lists is shown 
in the lower right. 
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Table 1 

Spitzer Extragalactic FLS MIPS-Ge Observations 



Field 




tint (s) 




5cr (mjy) 


Area 
(deg 2 ) 


70/im 


160/im 


70/im 


160/im 


Main 


42 


8 


19 


100 


4.5 


Verification . 


252 


42 


8 


50 


0.35 


Warm+Main 


84 




14 




3.3 



Note. — The integration time (t»„t) is the average integration time for the field. The approximate average 5a noise 
levels are for point sources and do not include the absolute flux density calibration uncertainties. 



Table 2 

Extragalactic FLS 70 \im Catalog 



Source Name 


a(J2000) 


<5(J2000) 


Error 


S70 


Flag 




(deg) 


(deg) 


(arcsec) 


(mjy) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


FLS70 J172723. 3+601626 


261.84712 


60.27414 


1.9 


81.4+12.9 


prf 


FLS70 J172716.2+600051 


261.81791 


60.01433 


1.6 


1610.1+244.6 


al20m 


FLS70 J172709.5+600141 


261.78960 


60.02815 


2.1 


77.3+12.7 


prf 


FLS70 J172704.7+593736 


261.76978 


59.62680 


2.0 


62.1+9.9 


prf 


FLS70 J172651. 1+601851 


261.71329 


60.31419 


3.3 


34.7+7.0 


prf 



Note. — The entire xFLS 70/im catalog (> la) is presented in the electronic edition of the Astronomical Journal. A 
portion of Table 2 is shown here for guidance regarding the form and content of the catalog. Column (1) gives the source 
name following the IAU designations (Lortet, Borde, & Ochscnbein 1994). Column (2) and column (3) are the right 
ascension (a) and declination (S) J2000.0 source positions in decimal degrees. Column (4) provides the radial positional 
uncertainty in arcsec (la). Column (5) gives the flux density measurement in mjy. The flux density errors include the 
uncertainties of the absolute flux density scale. Column (6) provides a flag for the flux density measurement method. A 
flag of "prf" indicates that the source was fitted by the PRF, while a flag starting with an "a" indicates that an aperture 
measurement was used (for cases of extended sources and/or close blends). The number following the "a" gives the 
aperture diameter in arcsec; e.g., al20 implies an aperture diameter of 120". Sources with aperture flags ending with a 
"m" represent bright galaxies not well fitted by the PRF and with nearby companions; in these cases, the flux density of 
the bright galaxy was derived by subtracting the PRF measurement (s) of the nearby faint point source(s) from the total 
aperture flux density. 
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Table 3 

EXTRAGALACTIC FLS 160 [im CATALOG 



Source Name 


a(J2000) 


<5(J2000) 


Error 


S160 


Flag 




(deg) 


(deg) 


(arcsec) 


(mJy) 




(1) 


(2) 


(3) 


' (4) 


(5) 


(6) 


FLS160 J172902. 1+600523 


262.25801 


60.08859 


3.4 


294.8±76.0 


prf 


FLS160 J172855.5+600026 


262.23039 


60.00619 


7.4 


154.4±44.2 


prf 


FLS160 J172851.3+600821 


262.21295 


60.13807 


3.8 


293.7±76.4 


prf 


FLS160 J172826.5+600539 


262.10963 


60.09320 


3.0 


1211.0±310.5 


al92 


FLS160 J172817.7+600100 


262.07271 


60.01573 


5.9 


188.4±50.97 


prf 



Note. — The entire xFLS 160/im catalog (> 7a) is presented in the electronic edition of the Astronomical Journal. A 
portion of Table 3 is shown here for guidance regarding the form and content of the catalog. Column headers are the 
same as described in Table 2. 



Table 4 

Extragalactic FLS 70/um Number Counts 



Average S„ 
(mJy) 


Slow 

(mJy) 


Shigh 

(mJy) 


Number 


Completeness 


dN/dSS 25 
(gal sr _1 Jy 1 ' 5 ) 


8.1 


7 


9 


46 


0.50±0.15 


2556±1030 


10.6 


9 


13 


61 


0.80±0.08 


2073± 575 


16.5 


13 


20 


40 


1.0±0.05 


1866± 522 


24.5 


20 


30 


21 


1.0±0.03 


1847± 582 


46.0 


30 


70 


18 


1.0±0.03 


1912± 626 


13.9 


12 


16 


230 


0.60±0.15 


2172±744 


18.4 


16 


22 


230 


0.90±0.08 


1883±453 


26.7 


22 


34 


200 


1.00±0.05 


1931±466 


40.3 


34 


50 


89 


1.00±0.03 


1800±451 


61.2 


50 


80 


54 


1.00±0.03 


1660±439 


110.9 


80 


140 


25 


1.00±0.03 


1698±514 


227.9 


140 


440 


17 


1.00±0.03 


1398±464 



Note. — The first five entries are the 70/Ltm counts for the verification field, while the remaining entries are for the 
warm+main 70/im survey. The uncertainties for the differential counts (dN/dS S 2 ' 5 ) include the Poisson noise, the 
completeness uncertainty, and the absolute flux density calibration uncertainty of 15%. The counts are based on the > 5er 
source lists, corrected for completeness. 



16 



Frayer et al. 



Table 5 

extragalactic fls 160^to number counts 



Average S„ 


Slow 


Shigh 


Number 


Completeness 


dN/dSS 2b 


(mJy) 


(mJy) 


(mJy) 






(gal sr -1 Jy 15 ) 


68.0 


56.4 


77.0 


17 


0.55±0.15 


17156± 8978 


85.9 


77.0 


97.5 


15 


0.85±0.08 


17541± 8155 


116.4 


97.5 


138.6 


13 


0.97±0.05 


14231± 6678 


101.8 


92.4 


112.9 


60 


0.55±0.15 


12682± 6104 


124.0 


112.9 


133.4 


45 


0.65±0.1 


13170± 5688 


152.8 


133.4 


174.5 


44 


0.80±0.08 


8825± 3674 


211.7 


174.5 


256.6 


43 


0.92±0.05 


8472± 3460 


320.7 


256.6 


420.8 


25 


1.0±0.05 


6399± 2738 


679.0 


420.8 


882.6 


10 


1.0±0.05 


5935± 2926 



Note. — The first three entries are the 160/xm counts for the verification field, while the remaining entries are for the 
main 160/xm survey. The uncertainties for the differential counts (dN/dS S 2 ' 5 ) include the Poisson noise, the completeness 
uncertainty, and the absolute flux density calibration uncertainty of 25%. The counts are based on the > 5a source lists, 
corrected for completeness. 



